Wave systems with direct processes and localized losses or gains: the non-unitary 
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We study the scattering of waves in systems with losses or gains simulated by imaginary potentials. 
This is done for a complex delta potential that corresponds to a spatially localized absorption or 
amplification. In the Argand plane the scattering matrix moves on a circle C centered on the real 
axis, but not at the origin, that is tangent to the unit circle. From the numerical simulations it 
is concluded that the distribution of the scattering matrix, when measured from the center of the 
circle C, agrees with the non-unitary Poisson kernel. This result is also obtained analytically by 
extending the analyticity condition, of unitary scattering matrices, to the non-unitary ones. We 
use this non-unitary Poisson kernel to obtain the distribution of non-unitary scattering matrices 
when measured from the origin of the Argand plane. The obtained marginal distributions have an 
excellent agreement with the numerical results. 
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I. INTRODUCTION 

In recent years there has been an intense study of sev- 
eral properties of classical wave chaotic systems, such as 
microwave cavities [H-IH and graphs elastic, acoustic, 
and optical resonators [5-8]. These systems inevitably 
exhibit energy dissipation (absorption) and many exper- 
imental and theoretical investigations have focused on 
the effect of this absorption on the scattering properties 
(see Refs. [1, and references therein). 

In addition to losses of energy the scattering properties 
are also affected by the imperfect coupling of the anten- 
nas to the system, that gives rise to a prompt response 
due to direct reflections [l2|. It has been shown that 
for chaotic systems the distribution of the sub-unitary 
scattering matrix, S, is modulated by a generalization of 
the Poisson kernel, Poisson's kernel squared, in a single 
port, or one channel, configuration (l3| . For an arbitrary 
number of channels, it was found that the Poisson kernel 
squared is the Jacobian of the transformation between a 
non-unitary scattering matrix with direct processes and 
a non- unitary one without such processes [14] . 

The Poisson kernel was developed first for unitary scat- 
tering matrices. In the one channel (one-dimensional) 
case, it is obtained from the analytical structure of the 
scattering matrix S and can be interpreted as the proba- 
bility to find S in a certain region of its space [l5l - [r7| . For 
an arbitrary number of channels the Poisson kernel is the 
probability density of S with maximum information en- 
tropy, where the only physically relevant parameter is the 
energy average S, known as the optical S matrix. In the 
optical model [l8[ S quantifies the direct processes fl7j ; 
in the absence of such processes, S — and the proba- 
bility density of S is a constant. 
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In a stationary random processes the analyticity- 
ergodicity conditions yield the Poisson kernel as the prob- 
ability distribution of an ensemble of S matrices de- 
termined by the ensemble average (S) [lTj- This en- 
semble represents an ensemble of macroscopically iden- 
tical systems that describes the statistical fluctuations 
of scattering pro perties in chaotic cavities with direct 
processes [l{§ [29- When (S) = 0, S is uniformly dis- 
tributed. Also, the Poisson kernel appears as the Ja- 
cobian of a transformation between scattering matrices, 
with and without direct processes [2~l[. A physical real- 
ization of that transformation was given in Ref . [22| ; this 
realization was helpful to demonstrate the equivalence, 
in certain limit, between the imaginary potential model 
and the voltage probe model, both used to implement 
dephasing or absorption in ballistic chaotic cavities [23| . 

Our purpose in this paper is to deep in the under- 
standing of the Poisson kernel. We study the non-unitary 
scattering matrix, S, in a one-dimensional (ID) problem 
with spatially localized losses or gains, in the presence 
of promt responses. We analyze the analytical structure 
of S to find its probability of stay in a certain region 
of its space. It is worth mentioning that local losses or 
gains have been of theoretical and experimental interest 
in chaotic systems. In this way, local losses have been 
used to model surface absorption in chaotic cavities [13] ; 
quasimodes in elastic cavities were measured when the 
absorption increase in a progresive way 25]. Also, lo- 
cal gains were used to select scar modes in multimode 
cavities (26[. Furthermore, a waveguide terminated by 
a perfect absorber, was used to study fidelity in chaotic 
microwave cavities [27|. By the other hand, the nearest 
level spacing statistics were investigated in open chaotic 
systems, as a function of the coupling to the environ- 
ment [28]. However, we recall that in the present work 
we leave chaos behind and we concentrate in the local 
character of the absorption or amplification. 

We organize the paper as follows. In the next section 
we use the analytical structure in the unitary case to 
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get the analyticity condition for non-unitary scattering 
matrices, in the one channel situation. This is done by 
adding an imaginary part to the energy. Also, we ap- 
ply that analyticity condition to non-unitary scattering 
matrices with constant modulus to obtain the Poisson 
kernel. In Sec. IIIII we consider a ID problem with lo- 
cal absorption or amplification and find the probability 
of stay of the scattering matrix S that describes it. We 
present our conclusions in Sec. IIVI 



II. 



ANALYTICITY CONDITION AND 
POISSON'S KERNEL 



The scattering matrix S(E) of systems where the flux 
is conserved is a unitary one that depends on the energy 
E. In the one channel case it is a I x 1 matrix that 
moves on the unit circle as E is varied, but it does not 
visit all points with the same frequency when the energy 
average S{E) of S(E), does not vanish. Using the fact 
that S(E) is an analytic function in the upper half of the 
complex- E plane, it is shown that the energy average of 
the m-th power of S(E) coincides with the m-th power 
of 03, 



S m (E) = S(E) 



(1) 



This result is known as analyticity condition. 

Dissipation or amplification can be modeled by adding 
an imaginary part to the energy |23l |29| , in which case the 
scattering matrix, denoted by S(E), is no longer unitary, 
it is sub-unitary or over-unitary, respectively. That is, 
S(E) can be obtained from an unitary S(E) by extending 
the energy to the complex plane, 



S(E) = S(E±iV ), 



(2) 



where Vb > and the plus (minus) sign holds for absorp- 
tion (amplification). Using Eqs. ((T|) and (J2J) it is easy to 
prove that 



S m (E) = S(E) 



(3) 



This means that, in the non- unitary case, the scattering 
matrix satisfies the same analyticity condition as in the 
unitary one. 

In the lxl case S is a complex number that can be 
written in a polar form as 



S = VRe l ° 



(4) 



where R is the reflection coefficient; R < 1 for absorption 
and R > 1 for amplification. R and 9 vary with E and, 
therefore, S moves in the Argand plane in the unit disk 
in presence of absorption, or outside of it for amplifica- 
tion. Following Ref. [T3| , we assume the existence of the 
measure 



dP(S) = p(R,6)dRd6, 



(5) 



which is the probability to find R between R and R + dR 
and 9 between 9 and 9 + d9. 

A very particular situation is the one in which the 
non-unitary scattering matrix moves along the circle with 
constant radius, as a function of E; we denote this lxl 
matrix by S' = \[R'e %e , where R! takes the fixed value 
Rq. Assuming that the energy average of S' is not null, 
we want to determine the probability to find 9' in the 
interval between 9' and 9' + d9' . In this case, 



with 



dP(S') =p'(R',0')dR'dO' 



p'(R',9') = 5(R> -R )q'(9'). 



(6) 



(7) 



This can be used to evaluate the average of the mth power 
of S' as follows, 



S' dP'(S') = < 



m/2 



e m0 q'(9')d9'. (8) 



Since q'(9') is a periodic function of 9' it can be expanded 
in a Fourier series, 
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(9) 



where the c m 's are constants that satisfy c m = c*_ m to en- 
sure q'{9') to be real. Using the analyticity condition ^ 
it is easy to see that 
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(10) 



With these coefficients the sum in Eq. ((9j> can be done to 
obtain 



Rq 



S" 



and, therefore, 

P '(R',9') 



2tt 



R' 



(11) 



5" 



2tt 



S'-S' 



■5(R'- Ro ) 



(12) 



Equation (|11[) is the same expression as that obtained for 
a disk of radius Rq in the problem of transfer of heat [3(| ■ 
Eq. ()12|) reduces to the Poisson's kernel of Ref [13] in the 
unitary case when we integrate it over the variable R for 
R = 1. 



III. LOCAL ABSORPTION (AMPLIFICATION) 
IN ID CAVITIES 

A simple model of a cavity with absorption or am- 
plification and direct reflection consists of a Dirac delta 
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FIG. 1. (a) In the presence of absorption the motion of S 
describes a circle of radius \/ r Ro < 1 displaced along the real 
axis with respect to the unitary case (v = 0). (b) The reflec- 
tion coefficient R as a function of ka takes values between a 
minimum and unity, (c) The phase 9' of the S' , seen from 
the center of the circle, shows the resonances of the system, 
(d) The distribution of 8' shows an excellent agreement with 
Poisson's kernel, Eq. <|11|) . 



function has a node just in the position of the delta po- 
tential. The real part of the potential affects only the 
distribution of the phase 9 along the circle. To illustrate 
this, we present in Fig. [T]the results for the absorption 
case with ua — va = 10 3 and ka varying from 10 4 up to 
obtain 35 resonances (the amplification case is quite sim- 
ilar). In Fig. [lja) we see that the motion of S describes 
a circle of radius v^o < 1- In the panel (b) of the same 
figure we observe that R varies between a minimum value 
Rmm = {2\/Ro — I) 2 and i? max = 1, where R « 0.91. 



If we translate S to S' = V R'e %0 such that R' becomes 
equal to the constant Rq, 9' behaves as shown in Fig. [3c). 
What it is interesting here is that the probability distri- 
bution of 9' is given by Poisson's kernel, as can be seen in 
Fig. [Ud), where we compare the numerical experiment 
with the theoretical result (ITTj) with S' taken from the 
experimental data. Therefore, two relevant parameters 

are needed, Rq and S'. As can be seen, the agreement is 
excellent. 

The probability distribution of S is easily obtained 
from the one for S' if we choose S as 



S = S' - ( 1 - 



(15) 



potential with complex intensity, located at a distance a 
in front of an impenetrable barrier. The potential is 



V(x) 



CO, 

(u =F iv) S(x 



x < 
a), x > 



(13) 



where u and v are positive constants; the minus (plus) 
sign corresponds to absorption (amplification). Notice 
the local character of the absorption or amplification. 
The incident waves to the potential suffer multiple scat- 
tering before they leave the cavity formed between the 
infinite barrier and the delta potential. The outgoing 
plane wave amplitude is related to the incoming one by 
the lxl scattering matrix given by 



S 



[sin ka + (k/a) cos ka] + i(k/a) sin ka 
[sin ka + (k/a) cos ka] — i(k/a) sinfca ' 



(14) 



where a = 2m(u ± iv)/H 2 and k = y / 2mE / h 2 . It is easy 
to see that S is a complex number which can be written 
in polar form as in Eq. (j4]), where 9 is twice the phase 
shift plus 7T and R — S^S, is the reflection coefficient; 
R < 1 for absorption and R > 1 for amplification. When 
v = 0, a = 2mu/h 2 and the unitary case is recovered 
(R = l). 

For the potential of Eq. ([T3|) . due to the imaginary 
part of the delta potential, the motion of S in the Ar- 
gand plane describes a circle of radius y/Ro ^ 1 displaced 
along the real axis. In this case R is not fixed but it is 
distributed in a certain interval. This circle touch the 
unitary one in the point R = 1 and 9 = tt where the 
delta potential is totally transparent because the wave 



with i?o the corresponding value for absorption or ampli- 
fication. From Eq. (TT2|) we get 



p(R, 



1 



R'{R,8) - 



S + (1 - s/W, 



271 \s-s 

x S[R'(R, 9) - R ] , 



(16) 



where the relevant parameters Rq and S are obtained 
from the experimental data. To get the marginal distri- 
butions, p(R, 9) of Eq. (fTB"]) can be integrated over one of 
the variables R and 9. Integration over 9 gives 



w(R) 



i?0 



s + (i - 



(l-R) 



R-(2^-l) 



R 



Re IS 



R 



x i ±9 



I + E^k)s + S 2 



(17) 



R 



T6(i?-1) , 



where 9(x) is the Heaviside function, Re(z) stands for 
the real part of z, and the upper (lower) sign corresponds 
to absorption (amplification). 

In similar way, integrating over R the marginal distri- 
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FIG. 2. The probability distribution of R, w(R), is shown 
on the left panels for absorption (upper) and amplification 
(lower). Note that the range of R is different for absorption 
than that of amplification, although the distribution has the 
same form. The distribution q(6) of 8 is shown on the right 
panels. It is the same for both cases. 



bution q(0) of 9 can be obtained: 
x R - S+(l 

Z7T 



^5) 



Ra-[l- V r R u ') 2 sin 2 i 



(18) 



where 



cost 



Ro-[l- 



sm 



The comparison of the marginal distributions w(R) 
and q(8) given in Eqs. fp~T|) and (fl8|). with the distri- 
butions coming from the numerical experiment for both, 
absorption and amplification, are shown in Fig. [2] for the 
same parameters of Fig.[TJ Again an excellent agreement 
is obtained. 



IV. CONCLUSIONS 

Extending the analyticity conditions to the non- 
unitary case, the non-unitary Poisson kernel was ob- 
tained. This was done by adding (substracting) an 
imaginary part to the energy for absorbing (amplifying) 
systems. As a physical realization, we studied a one- 
dimensional problem in which the lxl scattering ma- 
trix is non-unitary. The absorption (amplification) is lo- 
cated just at one position since the scattering potential 
was taken as a delta potential with complex intensity in 
front of an impenetrable barrier. In this one-dimensional 
cavity with local absorption (amplification) the S'-matrix 
moves on a circle with radius smaller (larger) than one; 
the center of the circle is not located at the origin of the 
complex plane, as happens in the unitary case. In a par- 
ticular reference framework, such as the one in which the 
module of the scattering matrix is constant, the phase of 
the S'-matrix is distributed according to the non-unitary 
Poisson kernel. The relevant parameters can be obtained 
from the numerical simulation and are: (a) the energy 
average of the scattering matrix and (b) the reflection 
coefficient. 

We are conscious that the model presented here is 
somehow artificial and, therefore, it presents some limi- 
tations. For example a delta potential is very difficult to 
realize in an experiment. Fortunately, one-dimensional 
elastic systems are good candidates to simulate one- 
dimensional quantum systems [3l| . In this sense the scat- 
tering matrix of a rod with a narrow notch, in which an 
absorbent foam is added as in Ref. [29 ]. could be closely 
distributed according to the non-unitary Poisson kernel. 
Finally, we expect that our one-dimensional model stimu- 
lates further investigations in two-dimensional problems, 
like chaotic cavities with local losses or gains. 
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